Measurement of the Total Active Solar Neutrino Flux at the Sudbury Neutrino Observatory with 

Enhanced Neutral Current Sensitivity 
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The Sudbury Neutrino Observatory (SNO) has precisely determined the total active (Vj.) *B solar neutrino flux 
without assumptions about the energy dependence of the Ve survival probability. The measurements were made 
with dissolved NaCl in the heavy water to enhance the sensitivity and signature for neutral-current interactions. 
The flux is found to be 5.21 ± 0.27 (stat) ± 0.38 (syst) x 10* cm^'s"' , in agreement with previous measurements 
and standard solar models. A global analysis of these and other solar and reactor neutrino results yields Am^ = 
T-1*qI X 10"' eV^ and 6 = 32.5^23 degrees. Maximal mixing is rejected at the equivalent of 5.4 standard 
deviations. 



PACS numbers: 26.65. +t, 14.60.Pq, 13.15.+g, 95.85.Ry 

The Sudbury Neutrino Observatory (SNO) Ot detects ^B 
solar neutrinos through the reactions 

V(, + c/ — > p -i- p -{- (CC), 
Vx + d -> p + n + v.v (NC), 
Vx + ^ v.v + e" (ES). 

Only electron neutrinos produce charged-current interactions 
(CC), while the neutral-current (NC) and elastic scattering 
(ES) reactions have sensitivity to non-electron flavors. The 
NC reaction measures the total flux of all active neutrino fla- 



vors above a threshold of 2.2 MeV 

SNO previously measured the NC rate by observing neu- 
tron captures on deuterons, and found that a Standard-Model 
description with an undistorted '^B neutrino spectrum and 
CC, ES, and NC rates due solely to Ve interactions was re- 
jected fz, '3"]. This Letter presents measurements of the CC, 
NC, and ES rates from SNO's dissolved salt phase. 

The addition of 2 tonnes of NaCl to the kilotonne of heavy 
water increased the neutron capture efflciency and the associ- 
ated Cherenkov light. The solution was thoroughly mixed and 
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a conductivity scan along the vertical axis showed the NaCl 
concentration to be uniform within 0.5%. 

The data presented here were recorded between July 26, 
2001 and October 10, 2002, totaling 254.2 live days. The 
number of raw triggers was 435,721,068 and the data set was 
reduced to 3055 events after data reduction similar to that 
in j|2|] and analysis selection requirements. Cherenkov event 
backgrounds from /3 - y decays were reduced with an effec- 
tive electron kinetic energy threshold T^ff > 5.5 MeV and a 
fiducial volume with radius Rfit < 550 cm. 

The neutron detection efficiency and response were cali- 
brated with a ^^^Cf neutron source. The neutron capture effi- 
ciency is shown in Fig. 1(a). The detection efficiency for NC 
reactions in the heavy water was 0.399 + 0.010 (calibration) ± 
0.009 (fiducial volume) for reif> 5.5 MeV and 7;fit<550 cm, 
an increase of approximately a factor of three from the pure 
D2O phase. Calibration of the detector's optical and energy 
response has been updated to include time variation of the wa- 
ter transparency measurements made at various wavelengths 
throughout the running period. A normalization for pho- 
ton detection efficiency based on ""N calibration data |4| and 
Monte Carlo calculations was used to set the absolute energy 
scale. '^N data taken throughout the running period verified 
the gain drift (approximately 2% per year) predicted by Monte 
Carlo calculations based on the optical measurements. The 
energy response for electrons is characterized by a Gaussian 
function with resolution cTj- = -0.145+0.392 VrJ+0.0353r^, 
where Te is the electron kinetic energy in MeV. The energy 
scale uncertainty is 1.1%. 

Neutron capture on ^^Cl typically produces multiple y rays 
while the CC and ES reactions produce single electrons. The 
greater isotropy of the Cherenkov light from neutron capture 
events relative to CC and ES events allows good statistical 
separation of the event types. This separation allows a pre- 
cise measurement of the NC flux to be made independent of 
assumptions about the CC and ES energy spectra. 

The degree of the Cherenkov light isotropy is reflected in 
the pattern of photomultiplier-tube (PMT) hits. Event isotropy 
was characterized by parameters /3i, the average value of the 
Legendre polynomial Pi of the cosine of the angle between 
PMT hits [5]. The combination /3i + 4^4 = /3i4 was selected 
as the measure of event isotropy to optimize the separation of 
NC and CC events. Systematic uncertainty on /5i4 distribu- 
tions generated by Monte Carlo for signal events was evalu- 
ated by comparing '^N calibration data to Monte Carlo cal- 
culations 1 6] for events throughout the fiducial volume and 
running period. The uncertainty on the mean value of /3]4 
is 0.87%. Comparisons of /3i4 distributions from '^N events 
and neutron events from ^^^Cf to Monte Carlo calculations 
are shown in Fig.^^b). The Monte Carlo calculations of /3]4 
have also been verified with 19. 8 -MeV y-ray events from a 
^H(p,y)*He source |7], high-energy electron events (domi- 
nated by CC and ES interactions) from the pure D2O phase 
of the experiment, neutron events following muons, and with 
low-energy calibration sources. 

Backgrounds are summarized in TableU Low levels of the 
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FIG. 1: (a) Neutron capture efficiency versus source radial position 
for the pure D2O phase (capture on D) and salt phase (capture on 
CI or D) deduced from a -'"Cf source, with fits to an analytic func- 
tion (salt) and to a neutron diffusion model (D2O). (b) Event isotropy 
from data and Monte Carlo calculations of a ^^^Cf source and an 
y-ray source. 

U and Th progeny ^''*Bi and ^oS'pj ^.^jj create free neutrons 
from deuteron photodisintegration and low-energy Cherenkov 
events from /3 - y decays. Ex-situ assays and in-situ analysis 
techniques were employed to determine the avenge Jevels of 
uranium and thorium during the experiment |3, 8, 9, 10]. Re- 
sults from these methods are consistent. For the ^^^^Th chain, 
the weighted mean of the ex-situ MnOx and the in-situ mea- 
surement was used. The ^^^U chain activity is dominated by 
Rn ingress, which is highly time-variable. The in-situ de- 
termination was used to estimate this background, because it 
provides the appropriate time weighting. The average rate of 
background neutron production from activity in the D2O re- 
gion is 0.72;^° 23 neutrons per day. 

Backgrounds from atmospheric neutrino interactions and 
^^^U fission were estimated with the aid of the code NU- 
ANCE 1 1 1] and from event multiplicities. 

^'^Na (which originates from neutron activation of ^■'Na) can 
also emit y rays which photodisintegrate the deuteron. Resid- 
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TABLE I: Background events. The internal neutron and y-ray back- 
grounds are constrained in tiie analysis. The external-source neutrons 
are reported from the fit. The last two Table entries are included in 
the systematic uncertainty estimates. 



TABLE IL Systematic uncertainties on fluxes for the spectral shape 
unconstrained analysis of the salt data set. f denotes CC vs NC anti- 
correlation. 



Source 


Events 


Deuteron photodisintegration 


7^ 1 +24.U 
'-'•^-23.5 


^ii(a, a)pn 


2.8 ± 0.7 


"■'**0(Q',n) 


1.4 ±0.9 


Fission, atmospheric v (NC -l- 




sub-Cherenkov threshold CC) 


23.0 ±7.2 


Terrestrial and reactor v''s 


2.3 ±0.8 


Neutrons from rock 


< 1 


^'^Na activation 


8.4 ±2.3 


n from CNO v's 


0.3 ±0.3 


Total internal neutron background 


111.3!^^^ 


Internal y (fission, atmospheric y) 


5.2 ± 1.3 


"'N decays 


< 2.5 (68% CL) 


External-source neutrons (from fit) 


84 5+^"-' 


Cherenkov events from y6 - 7 decays 


< 14.7 (68% CL) 


'AV events" 


< 5.4 (68% CL) 



ual activation after calibrations and ^'^Na production in the wa- 
ter circulation system and in the heavy water in the "neck" of 
the vessel were determined and are included in Table U SNO 
is slightly sensitive to solar CNO neutrinos from the electron- 
capture decay of '^O and '^F. 

Neutrons and y rays produced at the acrylic vessel and in 
the light water can propagate into the fiducial volume. Radon 
progeny deposited on the surfaces of the vessel during con- 
struction can initiate (a,n) reactions on '^C, '^O, and '^O, 
and external y rays can photodisintegrate deuterium. The en- 
hanced neutron capture efficiency of salt makes these external- 
source neutrons readily apparent, and an additional distribu- 
tion function was included in the analysis to extract this com- 
ponent (TablelUl. Previous measurements 12. 3^ .121 with pure 
D2O were less sensitive to this background source, and a pre- 
liminary evaluation indicates the number of these background 
events was within the systematic uncertainties reported. 

The backgrounds from Cherenkov events inside and outside 
the fiducial volume were estimated using calibration source 
data, measured activity, Monte Carlo calculations, and con- 
trolled injections of Rn into the detector. These backgrounds 
were nearly negligible above the analysis energy threshold 
and within the fiducial volume, and are included as an un- 
certainty on the flux measurements. 

A class of background events identified and removed from 
the analysis in the pure D2O phase ("AV events") reconstruct 
near the acrylic vessel and were characterized by a nearly 
isotropic light distribution. Analyses of the pure D2O and salt 
data sets limit this background to 5.4 events (68% CL) for the 
present data. 

To minimize the possibility of introducing biases, a blind 
analysis procedure was used. The data set used during the 
development of the analysis procedures and the definition of 
parameters excluded an unknown fraction (< 30%) of the final 
data set, included an unknown admixture of muon-following 
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Energy scale 


-3. 7, +3. 6 


" 1 .0,+ 1. 1 


±1.8 


Energy resolution 


± l.z 


±U. 1 


±U. J 


Energy non-linearity 


±0.0 


-0.0, +0. 1 


±0.0 


l^adial accuiacy 


-3.0 +3.5 


-2.6 +2.5 


-2.6 +2.9 


Vertex resolution 


±0.2 


±0.0 


±0.2 


Angular resolution 


±0.2 


±0.2 


±2.4 


Isotropy mean t 


-3.4,+3.1 


-3.4,+2.6 


-0.9,+l.l 


Isotropy resolution 


±0.6 


±0.4 


±0.2 


Radial energy bias 


-2.4,+1.9 


±0.7 


-1.3,+1.2 


Vertex Z accuracy t 


-0.2,+0.3 


±0.1 


±0.1 


Internal background neutrons 


-1.9,+1.8 


±0.0 


±0.0 


Internal background y's 


±0.1 


±0.1 


±0.0 


Neutron capture 


-2.5,+2.7 


±0.0 


±0.0 


Cherenkov backgrounds 


-1.1, +0.0 


-1.1, +0.0 


±0.0 


"AV events" 


-0.4,+0.0 


-0.4,+0.0 


±0.0 


Total experimental uncertainty 


-7.3,+7.2 


-4.6,+3.8 


-4.3,+4.5 


Cross section 113] 


±1.1 


±1.2 


±0.5 



neutron events, and included an unknown NC cross -section 
scaling factor. After fixing all analysis procedures and param- 
eters, the blindness constraints were removed. The analysis 
was then performed on the 'open' data set, statistically sepa- 
rating events into CC, NC, ES, and external- source neutrons 
using an extended maximum likelihood analysis based on the 
distributions of isotropy, cosine of the event direction rela- 
tive to the vector from the Sun (cos 9q), and radius within the 
detector. This analysis differs from the analysis of the pure 
D2O data 1^ 13 fl2il since the spectral distributions of the ES 
and CC events are not constrained to the shape, but are 
extracted from the data. The extended maximum likelihood 
analysis yielded 1339.6+^]j CC, n0.3^_lll ES, 1344.2^^^;^ 
NC, and 84.5^33 ^ external-source neutron events. The sys- 
tematic uncertainties on derived fluxes are shown in Table Hll 
The isotropy, cos Oq, and kinetic energy distributions for the 
selected events are shown in Fig.|2 with statistical uncertain- 
ties only. A complete spectral analysis including the treat- 
ment of differential systematic uncertainties will be presented 
in a future report. The volume-weighted radial distributions 
[p - (^fit/600 cm)-'] are shown in Fig.|3] 

The fitted numbers of events give the equivalent 
fluxes QIIl (in units of 10^ cm ^s '): 



1.59!°:°?(stat)!0:0^(syst) 



0.08V 

-o.26V^^"^^ ±0.10(syst) 
5.21 ±0.27 (stat) ±0.38(syst), 



2.21+°-?i(stat) 



and the ratio of the '^B flux measured with the CC and NC 
reactions is 



= 0.306 + 0.026 (stat) + 0.024 (syst). 



Adding the constraint of an undistorted ^B energy spectrum 
to the analysis yields, for comparison with earlier results (in 
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units of 10^ cm 



fir 



1.70 



:0.07(stat)!0:'»(syst) 



0.24 (stat):::°:^^(syst), 



2-13;^o;28(stat)_oQg 
4 9Q -1- n 1/1 ("otoA+O-ZS 



consistent with the previous SNO measurements J3 |^ 
[l^ . The difference between the CC resuhs above (with 
and without constraint on the spectral shape) is 0.11 
0.05 (Stat)!!? '"'' 



_o]5°(syst) X 10" cm 



-2s-l 



The shape-unconstrained fluxes presented here combined 
with day and night energy spectra from the pure D2O 
phase 1 12] place constraints on allowed neutrino flavor mix- 
ing parameters. Two-flavor active neutrino oscillation models 
predict the CC, NC and ES rates in SNO |17|. New anal- 
ysis elements included Super-Kamiokande (SK) zenith spec- 
tra 1 18|, updated Ga experiments 1 19j_20J, improved treatment 
of **B spectral systematic uncertainties, and CI and Ga cross- 
section correlations fiT]. This improved oscillation analysis 
applied to the pure D2O data reproduces other results l.21..22iil . 
A combined ;^'2 fit to SNO D2O and salt data i23ll alone yields 
the allowed regions in Am^ and tan^ 6 shown in Fig. 0] In 
a global analysis of all solar neutrino data, the allowed re- 
gions in parameter space shrink considerably and the LMA 
region is selected, as shown in Fig. Isja). A global analy- 
sis including the KamLAND reactor anti-neutrino results 1,2411 
shrinks the allowed region further, with a best-fit point of 
Am2 = 7.1+1-2 X IQ-^ eV^ and = 32.5!2-^ degrees, where 
the errors reflect 1 cr constraints on the 2-dimensional re- 
gion [Fig. |5lb)]. With the new SNO measurements the al- 
lowed LMA region is constrained to only the lower band at 
> 99% CL. The best-fit point with one dimensional projection 
of the uncertainties in the individual parameters (marginalized 
uncertainties) is Arn^ = ^.l+i'] x 10"^ eV^ and = 32.5![-^ 
degrees. This disfavors maximal mixing at a level equivalent 
to 5.4 cr. In our analyses, the ratio fe of the total flux to 
the SSM f25'l value was a free parameter, while the total hep 
flux was fixed at 9.3 x lO-' cm"^ s"'. 

In summary, we have precisely measured the total flux of 
active **B neutrinos from the Sun without assumptions about 
the energy dependence of the electron neutrino survival prob- 
ability. The flux is in agreement with standard solar model 
calculations. These results combined with global solar and re- 
actor neutrino results reject the hypothesis of maximal mixing 
at a confidence level equivalent to 5.4 cr. 

This research was supported by: Canada: NSERC, Industry 
Canada, NRC, Northern Ontario Heritage Fund, Inco, AECL, 
Ontario Power Generation, HPCVL, CFI; US: Dept. of En- 
ergy; UK: PPARC. We thank the SNO technical staff for their 
strong contributions. 
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FIG. 2: Distribution of (a) j8i4, (b) cos9q and (c) kinetic energy, 
for the selected events. The CC and ES spectra are extracted from 
the data using pu and cos 60 distributions in each energy bin. Also 
shown are the Monte Carlo predictions for CC, ES, NC + internal 
and external-source neutron events, all scaled to the fit results. The 
dashed lines represent the summed components. All distributions are 
for events with Tetf >5.5 MeV and Rflt< 550 cm. Diff'erential system- 
atic s are not shown. 
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FIG. 3: Volume-weighted radius (p) distribution (acrylic vessel=l). 
Shown are the distributions for CC, ES, NC+intemal and external- 
source neutrons, scaled from the fit, and separately determined exter- 
nal background distributions. 
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FIG. 4: SNO-only neutrino oscillation contours, including pure D2O 
day and night spectra, salt CC, NC, ES fluxes, with '*B flux free and 



hep flux fixed. The best-fit point is Am 
fe = 1.03, with xVd.o.f.=26.2/34. 
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day and night spectra, salt CC, NC, ES fluxes, SK, CI, Ga. The 
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^Vd.o.f.=70.2/8I. (b) Solar global + KamLAND. The best-fit point 
is Am- = 7.1 X 10"^ tan- 9 = 0.41, fg = 1.02. In both (a) and (b) the 
*B flux is free and the hep flux is fixed. 
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